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ABSTRACT 

Sterile neutrinos with masses in the keV range are viable candidates for the warm dark matter. 
We analyze existing data for the extragalactic diffuse X-ray background for signatures of ster¬ 
ile neutrino decay. The absence of detectable signal within current uncertainties of background 
measurements puts model-independent constraints on allowed values of sterile neutrino mass 
and mixing angle, which we present in this work. 


1 INTRODUCTION 

At present time there exists an extensive body of evidence that most 
of the matter in the Universe is composed of new, yet undiscovered 
particles - dark matter (DM). Observations of (i) galactic rotation 
curves, (ii) cosmic microwave background radiation, (iii) gravita¬ 
tional lensing, and (iv) X-ray emission of hot gas in galaxy clusters 
provide independent measurements of DM content of the Universe. 

Another major experimental discovery of the recent decade 
is that of neutrino oscillations. There are s eparate measu rements 
of neutrino oscillations in solar neut r inos lAhmad et alJ J2002)), 
atm ospheric neutr inos iFu kuda et all 1 1998th . and reactor neutri¬ 
nos tEguchi et alJ 120031) ). Neutrino oscillations can be explained 
if neutrino is a massive particle, contrary to the Standard Model 
assumption. This means that along with the usual left-handed (or 
active) neutrinos there may exist also right-handed or sterile neu¬ 
trinos. 

Conventional se a -saw mechanism _jMinkowskj Jl977j); 
Yanagidal tl979:l980ll:Fiell-Mann et alt tl 98(f):iRamond 1 1979k 
Mohamatr^&^enianovid tl98flt) : iGlashowt i\ 98(t) l of gener¬ 
ation of small active neutrino masses implies that the sterile 
neutrinos are heavy (usually of the order of GUT energy scale 
~ 10 10 - 10 15 GeV) and that their mixing with usual matter is of the 
order sin0 ~ 10~ 10 - 10 -15 . In addition of the smallness of neutrino 
masses, models of this type can explain baryon asym metry of th e 
Universe via thermal leptogenesis IFukugita & Yan agidal ll98filV ) 
and anomalous electroweak fermion number non-conservation 
fcuzmin et all {1 98 5 }) ). However they do not offer a DM candidate. 

Recently it was proposed that neutrino oscillations, the origin 
of the dark matter, and baryon asymmetry of the Universe can be 
consistently explained in the model called neutrino Minimal Stan- 
dard Model (vMSM) lAsaka et"ah . (: 2005 ): lAsaka & Shaposhnikovi 
(2005)). This model is a natural extension of minimal Standard 
Model (MSM), where three right-handed neutrinos are introduced 
into the MSM Lagrangian. In this extension neutrinos obtain Dirac 


masses via Yukawa coupling analogous to the other quarks and lep¬ 
tons of MSM and in addition Majorana mass terms are allowed for 
right-handed neutrinos. Unlike conventional see-saw scenarios, all 
of these Majorana masses (which are roughly equal to the masses 
of corresponding sterile neutrinos) are chosen such that the mass of 
the lightest sterile neutrino is in the keV range and the other ones 
are < 100 GeV — below electroweak symmetry breaking scale. In 
this model a role of the dark matter particle is played by the lightest 
sterile neutrino. 

The existence of a relatively light sterile neutrino has non¬ 
trivial observable consequ ences for cosmology and astrophysics. 
It was proposed in lDodelson & Widrow ( ll 993 ) that a sterile neu¬ 
trino with the mass in the keV range may be a viable “warm” DM 
candidate. The small mixing angle (sin# ~ 10~ 6 - 10 4 ) between 
sterile and active neutrino ensures that sterile neutrinos were never 
in thermal equilibrium in the early Universe and thus allows their 
abundance to be smaller than the equilibrium one. Moreover, a ster¬ 
ile neutrino with these parameters is important for the physics of 
supernova iFrver & Kusenkol 120051) ) a nd was proposed a s an ex ¬ 
planation of the p u lsar kick velocit ies iKus enko & Segrel 1 1 9971) : 
[Fuller et all f2003)); iBarkovich et a fl fcOM) ). 

In addition to the dominant decay mode into three active neu¬ 
trinos, the light (with mass m s < 1 MeV) sterile neutrino can decay 
into an active one and a photon with the energy E y = mj2. Thus, 
there exists a possibility of direct detection of neutrino decay emis¬ 
sion line from the s ources with big concentr ation of DM, e.g. from 
the galaxy clusters iAbazaiian et al. 1 2001 1:)). Similarly, the signal 
from radiative sterile neutrino decays accumulated over the history 
of the Universe could be seen as a feature in the diffuse extragalac¬ 
tic background light spectrum. This opens up a possibility to study 
the physics beyond the Standard Model using astrophysical obser¬ 
vations. 

Recently there has been a number of works devoted to 
the analysis of the possibility to disco ver sterile neutrino radia - 
tive decays from X-ray observations dAbazaiian et alJ 1 200l b): 
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Mapelli & Ferrara) {2005)). For example, it was argued by 
Abazaiian et alJ fcOOlbjt that if sterile neutrinos composed 100% of 
all the DM, one should be able to detect the DM decay line against 
the background of the X-ray emissi on from the Virgo cluster. Ac¬ 
cording to lAbazaiian et al. ‘ t (200lrf) the non-detection of the line 
puts an upper limit m s < 5 keV o n the neutrino mass (this limit was, 
however recently revised in lAbazaiianH 2005l). who finds th e restric ¬ 
tion m, < 8 keV) . It wa s also noted bv lAbazaiian et all fcOOlbl) : 
iManelli & FerrarJ <2005lf that one can obtain even stronger con¬ 
straints m s < 2 keV - from diffuse extragalactic X-ray background 
(XRB) under the assumption, that the dark matter in the Universe 
is uniformly distributed up to the distances, corresponding to_red 
shifts z <s 1, T ogether with the claim of Ha nsen et aj] J2002I) : 
fyieT et a mui), putting lower bound m s > 2 keV on the neutrino 
mass from Lyman a-forest observations, this would lead to a very 
narrow window of allowed sterile neutrino masses, if not exclude it 
completely. 

In this paper we re-analyze the limit imposed on the param¬ 
eters of sterile neutrino by the observations of the diffuse X-ray 
background (XRB). For that we are processing actual astrophysical 
data of HEAO-1 and XMM-Newton missions. There are several 
motivations for this. Namely 

(i) All the abov e restr iction s on sterile neu trino mas s 

jAbazajiat^Ud ] l2001bh: iDoleov & HanserJ <200211 : 

iMapelli & Ferraral <2005tl : I Abazai ianOOOsii ) are model de¬ 
pendent and based on the assumption that sterile neutrinos were 
absent in the early Universe at temperatures larger than few GeV. 
Depending on the model, the relation between the mass of sterile 
neutrino, the mixing angle and the present-day sterile neutrino 
density O s does change. In fact, to compute the sterile neutrino 
abundance one needs to know whether there is any substantial 
lepton asymmetry of the universe at the time of sterile neutrino 
production, what is the coupling of sterile neutrino to other parti¬ 
cles such as inflaton or super-symmetric particles, etc. 1 Moreover, 
even if these uncertainties were removed, the reliable computation 
of the relic abundance of sterile neutrinos happens to be very 
difficult as the peak of their production falls on the QCD epoch of 
the un iverse evolution, corr e spond ing to the temperature ~ 150 
MeV <Dodelson & Widrowl <1993h ). where neither quark-gluon 
nor hadronic description of the plasma is possible. Therefore, 
before the particle physics model is fully specified and the physics 
of hadronic plasma is fully understood, one can not put a robust 
restriction on one single parameter of the model such as m s . 

Therefore we aim in this paper at clear separation between the 
model independent predictions, based solely on astrophysical ob¬ 
servations and any statements that depend on a given model and 
underlying assumptions. To this end we treat m s and sin 9 as two 
independent parameters and present the limits in the form of an 
“exclusion plot” in the ( m s ;O s sin 2 26) parameter space. 

It should be stressed that our data analysis is not based on any 
specific model of sterile neutrinos and as such can be applied to any 
“warm” DM candidate particle which has a radiative decay chan¬ 
nel. In case of sterile neutrino the full decay width of this process 
is related to parameter sin 9 via Eq. J4} (see below). 

(ii) Contrary to the p revious works < Abazaiian et alJ <2001blf : 
iMapelli & Ferraral < 200A ), we argue that non-isotropy or “clumpy- 
ness” of the matter distribution in the nearby Universe does not 

1 For example, if the coupling of sterile neutrino to inflaton is large enough, 
the main production mechanism will be the creation of sterile neutrinos in 
inflaton oscillations rather than active-sterile neutrino transition. 


relax the limit on the neutrino mass. Indeed, the fact that signifi¬ 
cant part of the dark matter at red shifts z < 10 is concentrated in 
galaxies and clusters of galaxies just means that the strongest signal 
from the dark matter decay should come from the sum of the signals 
from the compact sources at z < 10. Taking into account that DM 
decay signal from z < 10 is some 2 orders of magnitude stronger 
than that of from z > 10, while the subtraction of resolved sources 
reduces the residual X-ray background maximum by a factor of 10, 
we argue that it would be wrong to subtract the contribution from 
the resolved sources from the XRB observations when looking for 
the DM decay signal. The form of XRB background spectrum with 
sources subtracted is, in fact, unknown and the assumption that it 
has a shape of initial spectrum, scaled down according to resolved 
fraction (as in Abazaiian et al .. i l2001blf ) requires additional justifi¬ 
cation. 

(iii) We find that more elaborate analysis of the data enables 
to put tighter limits on the allowed region of the parameter space 
(m s , O s sin 2 29) from the XRB observations. The idea is that the 
cosmological DM decay spectrum is characterized not only by the 
total flux but also by a characteristic shape. Being present in the 
XRB spectrum, it would produce a local feature with some clear 
maximum and a width greater than spectral resolution of the in¬ 
strument. Features of such a scale are clearly absent in the data 
(tinier features could be present in the spectrum due to e.g. element 
lines, but they can not produce a bump wider than spectral reso¬ 
lution). Therefore, one can find that adding to the standard broad 
continuum model of XRB the DM decay component in a wrong 
place results in decrease of the overall quality of the fit of the data 
by such a two-component model (increase of the^ 2 of the fit). The 
condition that the two-component model provides an acceptable fit 
to the data imposes an upper limit on the flux in the DM decay 
component which is much more restrictive than the limit following 
from the condition that the flux of the DM component should not 
exceed the flux in the continuum component. 

The paper is organized as follows. In Section 0 we compute 
the contribution of the radiative decay of sterile neutrino to the dif¬ 
fuse X-ray background (XRB) and compare its shape with that of 
measured XRB. We discuss the effects of non-uniformness of the 
DM distribution in Section rm In Section 0 we obtain a model- 
independent exclusion region from HEAO-1 and XMM-Newton 
observations. 


2 THE CONTRIBUTION OF DM DECAYS INTO THE 
XRB. 

The sterile neutrino with a decay width T and mass m s decays into 
an active neutrino and emits photons with a line-like spectrum at 
the energy E y = mj2. However, photons emitted at different cos¬ 
mological distances are redshifted on their way to the Earth so that 
as a result the photon spectrum is given by (see lMasso & Toldr3 
Il999t> : I Abazai ian et alJCOOIfll ): 

d 2 N Tn° DM 1 

dQ.dE 4n EH(m s /(2E) - 1) ' J 

Here n° DM is the DM number density at the present time and H{z) is 
the Hubble parameter as a function of red shift. The explicit form 
of H(z) depends on the cosmological parameters. This means that 
the expected dark matter decay spectrum is different for different 
cosmologies. We will be interested in z corresponding to the epoch 
more recent than radiation-dominated universe and, therefore, H(z) 
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Figure 1. Diffuse X-ray background spectrum from HEAO-1 mission . The 
black solid curve shows empirical fit Eq. 0 bv lGruber et al . i T99§). The 
reduced x 2 °f this fit is 1.2. Dashed (green and red) lines represent the 
result of fit of the same data to the model of the form (3 with added DM 
component. The DM decay component (dot-dashed curve) is calculated for 
the “concordance” model = §.l,£l matter = 0.3 and for m s = 36.5 keV. 
The green (long dashed) line represents a fit of DM with the mixing angle 
sin 2 26 = 1.9 • 10“ 12 and the best achievable fit has the reduced;^ 2 = 2. 
For the red (short-dashed) line the mixing angle is sin 2 2 6 = 2.4 • 10“ 13 and 
reduced^" 2 = 3, which we choose as a border line of allowed quality of fit. 


can be written as 

H{Z) * Ho V^A + 1 + z)3, (2) 

where Hq is the present-day value of the Hubble parameter and 
Qa , Quarter are the cosmological constant and matter contributions 
to the density of the Universe. Substituting □ into m we find that 
for small z the spectrum is approximated as 

<?N T n° DM (2 EV 2 ^ 

dQ.dE 2nH 0 yg£ 3 £2 A + Q mat , er m] ' 

Assuming the Majorana nature of neutrino mass, the decay width T 
is related to sin 20 via jPal & Wolfensteinl 1982j : iBareer et al.1199^1 : 

9aG\ , - n i / in, \ 5 , 

F = -- shr(20) m 5 = 5.6 x 10 ~ 2 sim 0 (--— ) sec 1 . 

256-4 n 4 s V1 keV I 

(4) 

As an example we show on Fig. Q the expected dark matter 
decay spectrum for the cosmological “concordance” model with 
Qa — 0.7, Qmatter — 0.3. One can see that close to the maximum the 
spectrum is roughly a power-law with the photon index equal to 1 

(dN/dE-E-f. _ _ 

The measurements of the diffuse XRB l Ma rshal et Ifl luicl) : 
[Gruber et al.Ul999i) l show that for 3 < E < 60 keV its form is well 
approximated by the following analytical expression 

d 2 F XRB , £ a E r r “ +1 keV 

- - C X rb exp(--—) — - 7 . (5) 

v T X rb ' 160 keV 1 keV ■ sr • sec • cm- 


In Fig.Qwe show the HEAO-1 data points 2 together with the above 
analytical fit (solid black line). Such a form of XRB spectrum can 
be explained by AG N emission, under c ertain as sumption s abo ut 
AGN populations jWorslev~ ID <2005t) ; iTreister & Urn! 1 20051). 
Below ~ 15 keV XRB was measured by many X-ray missions and 
the background was found to foll ow the powe r -law with the photon 
index V X m: - 1.3 - 1.4 jC^bere^lJjl^9jljLumbetan <20021) : 
iRevnivtsev etaD<20031 l2004 : lGillil<20031) i fBargeJJ2003I) ). Anal¬ 
ysis of iGrubeTetaTTriOoT ulso finds T xrb — 41.13 keV und 
Cxrb = 7.9. The DM decay component produces a harder spec- 
trum, as one can see from Fig.Q 

2.1 The uniformness of the DM density in the Universe 

Most of the power in the very hard DM decay spectrum of Fig.Qis 
emitted in the narrow energy interval close to the maximal energy 
E max — m s /2. From Eq. 0 one can see that emission in this en¬ 
ergy range is produced by neutrinos decaying at the present epoch 
(z - 0). For the case of the “concordance” model, the DM de¬ 
cay spectrum is characterized by the very hard inverted power law 
dN/dE ~ E l/2 below the energy ~ E max /2. The energy flux drops 
by an order of magnitude at the energies ~ E„ mx / 4 which corre¬ 
spond to the redshift z ^ 3. Thus, most of the DM decay signal is 
collected from the low redshifts. 

It is known that the process of structure formation leads to sig¬ 
nificant clustering of the dark matter at small redshifts. The signal 
from DM decays at small z is thus dominated by the sum of contri¬ 
butions from the point-like sources corresponding to the large DM 
concentrati ons, like g alaxies or cl usters of galaxies. 

lAbazaiian et all 1200 1 hi) and lMapelli & Ferran] 120051) argued 
that the clustering of the DM at small redshifts makes if difficult 
for the instruments that measure the diffuse X-ray background to 
detect the signal from the DM decays at small redshifts. Indeed, 
measurements of the diffuse X-ray background done with narrow- 
field instruments, like Chandra or XMM-Newton can miss the DM 
signal because of the absence of large nearby galaxies or galaxy 
clusters in the fields used for the deep observations and background 
measurements. As a measure for galaxy clustering one can take the 
distribution of the number of galaxies as a function of the distance 
(Afgai(r)). The Sloan Digital Sky Survey data shows that the function 
(Agai(r)) becomes constant for r > 1 00 Mpc (which corresponds 
to red shift z ~ 0.02) (see e.g. lJovce et alJ d2005t) ). Therefore, an 
instrument with the field of view (FoV), which encompasses the 
volume ~ (100 Mpc) 3 at distances, corresponding to z < 1, will 
observe homogeneous matter distribution. Let us take Hubble dis¬ 
tance as an estimate for such distances: Hf ~ 3.8 x 10 3 Mpc (here 
Hq is a present day Hubble constant). Then the minimal FoV 0f m , 
is determined from a simple relation 0 2 m Hf > (100 Mpc) 3 , i.e. 
ttfov <: 15'. The FoV of XMM-Newton is precisely of this order 
ttxmm ~ 30'. Therefore, the strongest signal from small z can be 
missed if only “empty fields” are selected for the XMM XRB ob¬ 
servations. 

However, wide field-of-view instruments or the instruments 
which have conducted all-sky surveys, like HEAO-1 or ROSAT 
cannot miss the largest contribution to the DM decay signal from 
z < 10 because of the full sky coverage. Thus, the above argument 
should be used in a reverse sense: to find the DM decay signal it is 
necessary to use the data on the X-ray background collected from 

2 We thank D. Gruber for sharing this data with us and for many useful 
comments. 
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m s [keV] 

Figure 2. Exclusion plot on parameters m s and sin 2 26 using HEAO-1 and 
XMM data. The values in the non-shaded region are allowed. In the re¬ 
gion where both HEAO-1 and XMM-Newton data is available, HEAO-1 
provides a more stringent constraint (as discussed in Section^^}. We sup¬ 
posed that sterile neutrino constitute 100% of all the DM. (i.e. Cl s = CIdm)- 
To remove dependence on the value of Qdm, we choose to plot £2, sin 2 26 
rather than sin 2 26. The results can be described by a simple empirical for¬ 
mula Eq.<6| (thin blue dashed line on the Figure). The dash-dotted green 
line represent exclusion region, obtained if one attributes all 100% of the 
XRB flux to DM at energies E <, m s /2. 

the whole sky, rather than front the “deep field” observations of a 
narrow-field instrument. 


3 RESTRICTIONS ON PARAMETERS OF STERILE 
NEUTRINOS FROM THE MEASUREMENTS OF THE 
DIFFUSE X-RAY BACKGROUND 

3.1 Restriction from HEAO-1 measurements 

The above analytical approximation J5} provides a good fit to the 
HEAO-1 data. No DM decay feature with the spectrum of the form 
of 0 (corrected for the spectral resolution of the instrument) is 
evident in the data. Straightforward constraint on the possible con¬ 
tribution of DM decays into the XRB spectrum is that the flux of 
the DM decay contrib ution does not exceed the total flux. Such an 
approach was used bv lDolgov & Hajisetj j20()A. who a pplied it to 
the broad-band XRB model of lResseHfeTumenll989t) . Using the 
fit j5j of lGruber et alJ (1999 ), which gives a much better descrip¬ 
tion in the keV region, we obtain an exclusion curve, shown as 
green dash-dotted line on Fig.0 

Clearly this restriction is model independent and rigorous. 
However, it can be made stronger. Indeed, the experimentally mea¬ 
sured XRB background is monotonic at energies from keV to GeV 
(see lOmber et al . 11999 )1. If the flux of the DM have composed 
significant part of all XRB, contributions to XRB, unrelated to 
the radiative DM decay, should have combined themselves into 
a spectrum with a sudden narrow “dip”. In particular, at ener¬ 
gies E > 111 J2 some new physical phenomena should have sud¬ 


denly “kicked in”, causing the spectrum to experience a very sharp 
“jump”. Moreover, the form of this dip must have been almost the 
same as the dark matter contribution with a minus sign. This would 
have required a mechanism of a very precise fine-tuning of con¬ 
tributions to XRB between DM and various physical phenomena 
or simply a chance coincidence. Of course, such a conspiracy can¬ 
not be absolutely excluded, in particular because of the lack of un¬ 
ambiguous theoretical predictions of the shape of XRB spectrum. 
However, it is very unlikely to have a precise cancellation of contri¬ 
butions of different nature. So, in what follows we will assume that 
such a situation does not happen (in particular, this does not hap- 
pen in models, attributing e xisting XRB shape to AGN emissions 
<Worslev et al] <2005ll : lTreister & Urrvl <20051) 1). 

Therefore, in the present work we argue that the constraint on 
the possible DM contribution to the XRB resulting from HEAO-1 
data can be improved using the statistical analysis of the data. Our 
strategy will be the following. We take the actual HEAO-1 data and 
fit it to the model of the form (varying C X rb , T X rb , f XRB ) plus 
an additional DM flux 0 (corrected for the spectral resolution of 
HEAO-1 instrument). Addition of a large DM contribution would 
worsen the quality of the model fit to the data, while addition of 
a small DM contribution does not change the overall fit quality. 
Thus, one can put a restriction on the DM contribution by allowing 
the DM component to worsen the fit by a certain value. Taking into 
account that we have around 40 degrees of freedom of the system 
under consideration, we take the maximal allowed value of the re¬ 
duced chi square of the fit to be x 2 < 3. Thus on technical level 
our method restricts the flux of DM to be of the order of errors of 
measured XRB flux, rather than its total value. 

The results of the application of the above algorithm to the 
data are shown on the Fig.0 The allowed values of (m s , £2, sin 2 26) 
are those in the unshaded region. 3 The shaded region below the 
dash-dotted line is excluded under the assumption that XRB spec¬ 
trum does not have an unknown feature (a “dip”), fine-tuned to be 
located precisely at energies where DM contributes. 

According to Eqs. 0 J4} flux of DM is proportional to 
£2 S sin 2 26. Therefore we choose to plot this value rather than con¬ 
ventional sin 2 26 on y axis of Fig.Q This removes uncertainty, re¬ 
lated to the determination of Q-dm- Notice, that the results obtained 
from HEAO-1 data (thick solid lines on Fig. [3 are model indepen¬ 
dent and can be applied to any DM candidate, which has a decay 
channel into a lighter particle and a photon, with its decay width T 
related to sin 9 via Eq. o 


3.2 Restriction from XMM background measurements 

In the energy band below 10 keV the XRB was studied by 
numerous narrow FoV instruments, including XMM-Newton. 
Better angular resolution and sensitivity of these instruments 
has ena bled to resolve some 90% of the XRB into poin t 
sources < Brandt & HasingeJ 1200-51) : iDe Luca & Moleiidil <2005l) : 
iBauer etan ~ l2004ll ). Apparently, better sensitivity should enable 
to put tighter constraints on the possible DM contribution into the 
XRB. However, in this section we show that this in not the case. 

3 The normalization of XRB spectrum, measuredby HEAO-1, is known 
to be lower than any other XRB measurements <Gilli 12003 1); iMoretti et all 
120031) 1. Therefore on F ig.[2lHEAO -l data were increase d by 40% accord¬ 
ing to lPe Luca & M olendFl20051): IWorslev et all l2005t) : iTreister & Urrvl 
1 2005) (see, however iGilli l2QQ3h ). This weakens the restriction of Fig. |2] 
by about 10%, as compared to actual HEAO-1 data. 
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The key point which leads to such a conclusion is that the bet¬ 
ter sensitivity (to point sources) of these instruments is achieved 
due to the better angular resolution, rather than due to the increase 
of the effective collection area. Thus, if one is interested in the dif¬ 
fuse sources, the sensitivity is not improved compared to HEAO-1. 
On the contrary, in the narrow FoV instruments the XRB signal is 
collected from a smaller portion of the sky, which leads to lower 
statistics of the signal. At the same time, the instrumental back¬ 
ground (which is thought to be due to the cosmic rays hitting the 
instrument) is roughly the same for narrow and wide field instru¬ 
ments. Thus, the ratio of the counts due to the photons of the XRB 
to the instrumental background counts is smaller for the narrow 
field instruments and one needs larger integration times to achieve 
good statistical significance of the XRB signal. 

The above problem would not have affected the constraints 
on m„ sin 6 if they were imposed by the condition for the flux in 
the DM decay component not to exceed the total XRB flux in a 
given energy interval. However, as the statistical significance of the 
XRB signal in the narrow FoV instruments is lower, bigger errors 
decrease x 2 value, and therefore the imposed limit on DM will be 
worse than that of the corresponding wide FoV instruments. 

The exclusion plot in the ( m s , Q s sin 2 28) parameter space ob¬ 
tained from the analysis of the XMM data, is shown in Fig.0 We 
took the a ctual data of two collections of XMM background ob- 
servationsjLumbeytL! (2002), total exposure time ~ 450 ksec and 
iRead & Ponmanld200 J) total observation time ~ 1 Msec). The form 
of the XMM background is fitted by the power-law with the index 
F = 1.4 jLumb etalJ<2002lP . in agreement with Eq. Applying 
to XMM data our method, as described in the previous Section, we 
find the restrictions, shown on Fig. [2] In the neutrino mass region 
6 keV < m s < 20 keV where both HEAO-1 and XMM data are 
available, the XMM-Newton background provides weaker restric¬ 
tion that the data from HEAO-1, because of the reason discussed 
above (the XMM FoV has the radius 15’, which is much smaller 
than that of the HEAO-1, ~ 3° for A2 HED detectors). In case 
of XMM the statistical error at energies above roughly 7 keV is 
dominated by instrumental background, while in case of HEAO-1, 
the errors at these energies were dominated by diffuse X-ray back¬ 
ground. As a result, errors of flux determination are smaller for 
HEAO data, thus providing more stringent restrictions. Of course, 
the XMM data provides the constraints on parameters of the neu¬ 
trino in the region 1 <, m, < 6 keV where HEAO-1 data are not 
available. The data on Fig.^fits to the simple empirical formula: 

n jS in 2 (20)<3xl(T 5 (^) 5 . (6) 


3.3 Discussion. 

In this work we have looked for signatures of sterile neutrino decay 
in the extragalactic X-ray background. The main result of our pa¬ 
per is the plot on Fig. |2| which constraints the properties of sterile 
neutrino as a dark matter candidate. 

Let us compare our exclusion plot with the similar plots, found 
in other works. |Solgg^^Hl^isgr| used a broad-band limit, 

put on XRB by lResse^j^umen Il98^) . to find a restrictions on 
sterile neutrino parameters. First of all, in the keV energy band the 
limit o f lRessel & Turned 119891) is weaker than that of Grub er et alJ 
(1999), which we are using in this paper. Therefore the dashed- 
dotted line on the Fig. 0 provides model-independent restrictions 
on parameters of sterile neutrino that are str onger than the similar 
bound on Fig. 4 in iDoleov & Hansenl 12002 !). Second, as discussed 


in Section[3] we put a restriction on parameters of sterile neutrino, 
based on the statistical analysis of the actual experimental data of 
HEAO-1 and XMM-Newton missions. Such a constraint is possi¬ 
ble under the assumption that there is no extremely unlikely fine- 
tuning between various components, contributing to the XRB. To 
put it differently, the observed XRB spectrum is monotonic in keV 
range (at energies below ~ 15 keV it was checked by various X-ray 
missions, including recent measurements by XMM and Chandra). 
Therefore, we assume that XRB spectrum without DM has no “dip” 
with the shape and location precisely at the place, where DM con¬ 
tributes. The results of our analysis are shown in thick black solid 
line on Fig.[^] As the statistical quality of the data is good, this puts 
about two orders of magnitude stronger restrictions then those, rep¬ 
resented by the dashed-dotted line. Region between the solid (or its 
empirical fit, Eq. JfJ) and dashed-dotted lines is excluded under the 
assumptions discussed above. 

Ref.lAbazaiian et~al ] feOOlbj ) provides an exclusion plot (also 
quoted by A bazajjarj j200fj )) base d on observations of the Virgo 
cluster. Ref lBowj^kveTan J2006ll shows that t he restrictions from 
clusters of galaxies are weaker than those of Abaz aiianl 12005I) . 
however, in mass ranges 2 keV < m s < 10 keV they are 2 to 4 
times better than those, coming from XRB. For careful comparison 
of our results with those of lAbazaiianl 120051) see lBovarskv et all 
fcOOfil) . 

I n some cases, such as in the Dodelson-Widrow (DW) sce¬ 
nario foodelson & Widrowl i 1 99 it) ) 4 there exists a relation be¬ 
tween the mass of sterile neutrino, the mixing angle and the 
present-day ster i le neutr i no density Q., {Dol gov & Hansenl 120021) : 
lAbazaiian et alJ l2001hl) : lAbazaiian J2005to . Since both the sterile 
neutrino abundance and the gamma-ray flux depend on the mixing 
angle 8 and sterile neutrino mass m s only, this provides an upper 
limit on sterile neutrino mass in the DW scenario. Accordin g to 
the computation of sterile neutrino abundance made in lAbazaiianl 
1 2005 ). the corresponding limit following from the X-ray bound 
found in the present paper reads m s < 8.9 keV for the DW ster¬ 
ile neutrino. However, even in DW scenario this number is subject 
to uncertainties which are difficult to estimate, since the peak of 
the sterile neutrino production falls into the range of temperatures 
where neither hadronic nor quark description of the hot plasma is 
possible. 

The Standard Model augmented by just one sterile neutrino 
cannot explain the neutrino oscillations data and thus extra ingre- 
dients should be added to it. The analysis of lAsaka et al. I 11 2006 ) 
shows that in a more realistic model, the vMSM, a relation between 
the mass of sterile neutrino, the mixing angle and the present-day 
sterile neutrino density £2 S changes, depending on unknown param¬ 
eters, such as the masses and couplings of extra particles and pri¬ 
mordial abundances of sterile neutrinos. Therefore, for a general 
case one cannot put a model-independent restriction on the mass of 
sterile neutrino m s , and the search for a narrow line coming from 
its decays should be done in all energy ranges. 

The results of this work show that a possible strategy to look 
for a feature of DM decay in the XRB is to use data of “all sky sur¬ 
veys” rather than “deep field observations” and utilize instruments 
with the largest possible FoV (rather than best angular/spectral res- 


4 In this scenario it is assumed that there is just one species of sterile neu¬ 
trinos, that their concentration was zero at temperatures higher than 1 GeV, 
that there were no substantial lepton asymmetries, and that the universe re¬ 
heating temperature was larger than 1 GeV. 
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olution). An example of such an instrument can be INTEGRAL, 
which currently plans an observation of XRB. 
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